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The Schrodinger equation is used to exactly evaluate the propagator, wave 
function, energy expectation values, uncertainty values and coherent state for a 
harmonic oscillator with a time -dependent frequency and an external driving time- 
dependent force. These quantities represent the solution of the classical equation 
of motion for the time -dependent harmonic oscillator. 

I. Introduction 

It is well known that an exact solution of the Schrodinger equation is 
possible only for special cases. For this reason, approximate methods are needed. 
Exact solutions provide important tests for these approximate methods and for 
various models of physical phenomena. In general, the solution of the Schrodinger 
equation for explicit time - dependent systems has met with limited success because 
of analytical difficulties, although progress has been made during the past three 
decades. 1 " 5 Camiz et al® have obtained the wave functions of a t ime - dependent 
harmonic oscillator perturbed by an inverse quadratic potential, using the 
Schrodinger formalism and a generating function. Further, Khandekar and Lawande 
have evaluated the exact propagator and wave function for a time -dependent harmonic 
oscillator, both with and without an inverse quadratic potential, using Feynman 
path integrals. In addition, Jannussis et al® have calculated the propagator for 
several quantum mechanical systems with friction. 

In a previous paper, ^ we have evaluated the propagator, wave function, energy 

expectation values, uncertainty values and transition amplitudes for a quantum 

damped driven harmonic oscillator by using path integral methods. Also, we have 

10 

obtained the coherent state for the damped harmonic oscillator and calculated the 
propagator for coupled driven harmonic oscillators. 
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In this paper we discuss the exact quantum theory of a forced harmonic 
oscillator with a time -dependent frequency. In Sec. II we evaluate the propagator 
using the Schrddinger equation and path integral methods, and in Sec. Ill we 
calculate the wave functions using the propagator. In Sec. IV we define the energy 
operator and calculate energy expectation values. In Sec. V we obtain the 
uncertainty values. In Sec. IV we determine the coherent state and its properties. 
Finally, In Sec. VII we present results and a discussion. 


II. Propagator 

We consider a system whose classical Hamiltonian is of the form 

H ' 2M p2 + 2 “ 2(t) 5(2 ‘ f(t)x • (2.1) 

where x is a canonical coordinate, p is its conjugate momentum, o(t) is a frequency 
as a function of time, M is a positive real mass, and f(t) is an external driving 
force. The Lagrangian corresponding to the Hamiltonian (2.1) is 

1*2 1 2 2 

L — 2 " 2 (t)x + f(t)x . (2 2) 

Here, the Hamiltonian H and Lagrangian L depend on time. The classical equation of 
motion for our system is 


d 2 1 

2 x + £J ( t ) x - - f(t) 


For the case where u(t) - (constant), the solution of Eq . (2.3) represents 
harmonic motion; otherwise, it is difficult to evaluate the exact solution. 

The path integral formulation of Feynman provides an alternate approach to 
solving dynamical problems in quantum mechanics. 12 In this approach, the usual 
Schrodinger equation is replaced by the integral equation 


V>(x,t) 


- J dx* K(x,t;x't') tf(x’.t') (t > f) 


with the initial condition *(x,t) - iKx'.t). Here, *(x,t) is a wave function and 
K(x,t; x'.t') is a propagator. The propagator K(x,t; x'.f) is defined by the path 
integral 


f(x,t) N-l 

K(x,t; x',t') - lim II dx exp[^ S(x,t; x',t') 

J (x' ,f) j-1 J »* 
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where the integration is over all possible paths from the point (x',t') to the 
point (x,t), and S(x,t; x',t') Is the action defined as 

S(x,t; x',t') - f dr L(x,x,r) (2.6) 

J t' 


For a short time interval £, substitution of Eqs . (2.2) and (2.6) into Eq . (2.5) 
gives the normalizing factor A and the usual Schrddinger equation: 


A. 

J 




(2.7) 


itf 


at 


*“5 * + \ Mu, 2 (t) X 2 V - f(t) x i> 
2M 3x 


( 2 . 8 ) 


Since K(x,t; x't') can be thought of as a function of the variables (x,t) or of 
(x',t'), it is a special wave function, and it satisfies Eq. (2.8): 

2 2 

IK ~ K(x.t; x',t') - - ^ ^ K(x,t;x'.t') + \ M^ 2 (t) x 2 K(x,t; x\t') 

dx 


- f(t) x K(x , t ; x f , t ' ) , (t > t') 


i\i — , K(x, t;x' , t ' ) - 


2 2 

^ — 2 K ( x , t ; x't') + ^ x ' 2 

dx' 


(2.9) 


x K(x , t ; x' , t ' ) - f(t') x ' K(x,t;x',t') 


(t' > t) . (2.10) 

* 12 , 13 

Because the Lagrangian is quadratic, the propagator has the rorm 

o 2 

K(x, t; x't') - exp[a(t,t')x + b(t,t , )xx < + c(t,t')x' + g(t,t')x 

+ h(t,t')x' + d(t,t')] , (2.11) 

where from Eqs. (2.9) and (2.10) we can easily deduce that the coefficient of the 
third and higher powers in x is zero. 

Substituting Eq . (2.11) into Eqs. (2.9) and (2.10), we obtain the differential 
equations 
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d 2iH 2 H 2. . 

— a - „ a + w (t) 

dt M 2 ip 


( 2 . 12 ) 


^ (bx'+g) - ^ a(bx'+g) + J f(t) 


(2.13) 


— (cx' 2 +hx'+d) " ^ a + 2M (bx ' +B)< 


(2.14) 


d 2M 2 m . t .v2 

C ‘ Mi C + 2* “ (t } 


(2.15) 


i (bx+h) - jjf c (bx+h) + 


(2.16) 


jh (ax 2 +gx + d) - ^ (bx+h) 2 + Jj c . 


(2.17) 


Equations (2.12) and (2.15) are nonlinear equacicr.c. For the cr.cc vhr-re w(t) - w q , 
a solution is easily found, but in other cases it is difficult to find an exact 
solution. If q(t) obeys the differential equation 


q(t) + c 2 (t) q(t) - 0 , 

dt 


(2.18) 


then the solutions of Eqs. (2 . 12) - (2 . 14) are 


iM q ( t) 

2tf q ( t) 


b ( t )x ' + g(t) 


ii,r 

V q L J 


ds f(s)q(s) + b 


(2.19) 


( 2 . 20 ) 


2 - H b Q r t ds b o f ds _ 

c(t)x' + h( t)x' + d(t) - Inq + 21M y J 2 iMM J , ,2 


dp f(p)q(p) + 


•■>*» j q(s) 2 j 


dp f(p)q(p) dr f(r)q(r) + d Q , (2.21) 


where b and d are constants of integration and do not depend on t, and the 
o o 

solutions of Eqs. (2 . 15 ) - (2 . 17 ) are 


c(t') 


M a(t f ) 
2iM q(t') 


( 2 . 22 ) 


b ( t ' ) x + h( t ' 


ds f (s)q(s) + b ' 
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• 2 ft' 


2 -h 1D o r c ds r t d S 

a(t')x + g(t')x + d(t') - Inq(t') + jjflT J 2 + K* J q(g) 2 

x | dp f(p)q(p) + | J d P f (P)^(P) J dr + d ; • (2-24) 


q(s) 


where b' and d' are constants of Integration and independent of t' . Since only t 
o o 

is a variable inEqs. (2 . 19) - (2 . 21) , we have suppressed t' in a(t,t'), b(t.t'), 
etc., and we have similarly suppressed t in Eqs. (2 . 22) - (2 . 24) . 

In polar form we may write 


q(t) - Jj(t)e 


i-r(t) 


(2.25) 


where rj(t) and 7 (t) are real quantities. From Eqs. (2.18) and (2.25) we note that 


fj(t) - r; ( t) 7 ^(t) + u 2 (t) rj(t ) - 0 
2»? (t)7(t) + fj(t) 7*(t) - 0 
'«? 2 (t)}(t) - 0 


(2.26) 

(2.27) 

(2.27') 


where the constant n is a time-variant quantity. From Eqs. (2.26) and (2.27), we 
find another form for the solution of Eq . (2.18) as 


q ( t) - »? ( t) sin( 7 - 7 # ) 
q ( t ' ) - r? ( t ' ) sin(7-7' ) 


(2.28) 

(2.29) 


where 7 - 7 (t) and 7 ' - 7 (t')* 

Substitution of Eq. (2.28) into Eqs. (2.19) and (2.21) gives 


a(t) " 9 5 i? + i cot(7-7'): 
rj 


»<«*■ * *"> - * nlU r * ,(s>f(s) sin[i ' <s) - 7 ' 1 


(2.30) 


(2.31) 


c(t)x ,Z + h(t)x' + d(t) - Intij* 1 * sin 4 (7-7')] + 2 ynM cot ('*‘T'') 


lb 
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iJttlM sin (7 


- r 


ds f(s)r?(s) sin[7(s) -7' ) ] 


2iHnM sln (7 


in(7-7' ) } dS | dp f ( s )f(p)»?(s)»»(p) 


X sin[7(s)-7' ]sin[7(p)-7') 


Furthermore, substitution of Eq . (2.29) into Eqs . (2.22) and (2.24) gives 


c(t') " 2)5 [ ‘ V + cot(T-'r' ) 1 


b(t')x + h(t') - 1 ^. sin ( 7 , 7 .) + l^rjsinC-r 


i r 

n( 7 " 7 # ) J 


ds f(s)r/(s)sin[7-7(s)] 


9 -h -U ib o 

a(t')x + g(t')x + d(t') - ln[r?' sin ( 7 - 7 ')] + cot ^ 7 ’ 7 ^ 


]rfOMsin(7-7' 


ds f(s)rj(s) sin[ 7*7(s) ] 


t # t 

n( 7 . 7 .) } ds | d P f(s)f(p)f?( 


2inHMsin(7 


s)n(v) 


X sin[7(s)-7]sin(7(p)-7] 


From Eqs. (2.31), (2.32), (2.34) and (2.35), we deduce that the constants 

and b^ are given as 


b - - M 7 ' f 7 ' x ' 
o 


b' - M7nx 
o 


Also, from the normalization condition, 


d o - ln ^ )H 


From Eqs. (2 . 30) - (2 . 37) and (2. 27'), we find that 
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a(t.t') - [J + 7 cot( 7 -7')] 


(2.39) 


“ 2 ^ + 7 * cot <7-7')] 


b(t t<) _ m -Jkxi — 

1 ' ; y sin(7-7') 

- K h £rr> C ds ' <■>-’' 

- » sin?v-y) £. ds ^ 

— — f ds -=^= sin[ 7 - 7 (s)] [ dp ~= 

,in(7 7 > Jf Ms) J t- y; 


(2.40) 


(2.41) 


(2.42) 


(2.43) 


d(t.t') - - 


2M>( si 


77(p) 


x sin[7(p)-7' 


(2.44) 


Inserting Eqs . ( 2 . 39 ) - ( 2 . 44) in Eq. (2.11), we obtain the propagator for the forced 
time - dependent harmonic oscillator as 


Kfx t . x . t .' l _ r — m LkJn — 

M ’ ’ ; |>1* sln( 7 - 7 ')_ 


, im ,n 2 ai , 2 . , 

x exp[- (^ x - Jr )] 


x exp { 2 H sin( 7 - 7 ') [ (7x2 + ^' x ' 2) cos (.1'l ' ) - 2777' xx' 

.2 ii f* 

M V y;<«) 

♦ 4 - f . 

J t' 


ds sln[ 7 (s)- 7 *] 


ds - sln[ 7 - 7 (s)] 


77(S) 


. i- f ds =£^ sin[ 7 -7(s)] f dp sin[ 7 (p)-7' ]]} , 

M J t' 7 * 7 (s) J t' 77(P) JJ 


(2.45) 


where the unprimed and the primed variables denote the quantities which are 

functions of time t and t ' , respectively. It may be easily verified that for the 

case where w(t) is a real positive constant , we have rj ( t) - 1 and 7 (t) - u t, 

r o o 
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and the propagator of Eq . ( 2 . 45 ) reduces to the usual expression for a forced 

12 

harmonic oscillator. 


III. Wave function 

We now rewrite the propagator in another form in order to derive the wave 
function : 


r • « K 

M (y'y) 

Xt) - [ 2*1K .ln( 7 - T ') , 


exp [3 {[v x2 + & x I ds i 

^ x ' 2 + jjyv x' J 


Usl 


Jy(s) 


c °s I 7 ' 7 (s) ] 


] 




ds -=£=1 cos[7' -7(s) )j|j 


Jy(s) 


6X P [2^ cot (y"r' )| 77 X - J | ds sin t 7 - 7 (s) ] j ‘ 

\jy' X* - jjj [ ds -Ull sin[7' -7(s) ] 1 2 | 

L J J 7 (s) J ^ 

M* sin(i-7') {[ y; X ‘ M 1 dS S ln l7-7 (s) ] 

v7( s ) 

7W' - if-- 15 * 1 


x exp 


L2KM 


ds -=±=± sin[7' -7(s) ] 1}1 
|cot(7-7') IT ds - — sin[7-7(s) ] 1 ' 

1 u y;( S ) J 


+ cot(7 




dp sin[ 7 '- 7 (p)]l 2 

/ftp) J 


sin(7 


- 7 ') I dS sin( 7 - 7 S ) J < 


ds sin(7-7 s ) | dp = ^ sin[7’-7(p)] 


Jl(p) 


J f 

( 7 - 7 ‘) J t , 


sl n (7-7-; j t , J-f(s) 

x sin[7(p) -7 ' ] || 


ds - sin[7-7(s)] 


J c- y 7 ( p ) 


( 3 . 1 ) 
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H,h , • , , k , • k e ' 1( '-’'I 


- (Jjj> W) <7) 


1 - e 


■21(7*7') 


eXP [Mfe x2 + M ^ X | ds -MMX cos[ 7 -7(s)]] 

V7(s) J 

X ' 2 + M | ds 7^ cos[7'-7(s)]1}1 

J J 7(s) J '* -* 

M [ ds 7“~ sin (7-7(s) ]1 2 

M J y;<s) J 



ds sin[ 7 '- 7 (s)] 12 

77(s) 


exp fl* x ■ « f ds ■ iniT - 7<,,i ] : 

y-y'x' - “ f ■*- — ^ 


ds sln[ ( 7 - 7 (s) ]1 ' 

77(8) J 


x ’ M f ds 7“ sin [7-7(s) ]1 |y 7 'x' 

MJ 74(s) JL 


-if 


£ ( s } 

ds sin[ 7 '- 7 (s)] 


77 (s) 


}] 


e -itf(t) e ifl(t*> 


where 


6(t') - 6(t) - 2 ^ |cot(7-7')|J ds sln t7-7(s) ]J‘ 


+ cot (7- 7 ' ) 


t' 


ds sin[ 7 f - 7 (s)] ^ + 2 ‘ ~ 


74(s) J ~ sin( 7 - 7 ') 

t' 


ft pt' 

x ds 7- — sin[ 7 - 7 (s) ] dp sin[ 7 '- 7 (p)] 


Jl (s) 


v7(p) 

- |" ds sin[ 7 - 7 (s) ] [ dp — ■ sin[ 7 (p)- 7 '] 

J t' J-r(s) J t' 74(p) 


Let us introduce Mehler's formula, 


14 


(3.2) 


(3.3) 
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(3.4) 


exp [ - (X + Y 2 - 2XY) / ( 1 - Z } 1 g -(xW) X „.. Z n . H (y) 

Jl - Z 2 J_ Q 2 n n! n n 


where 


X - [77 x + jjj f ds -Mil sin[7-7(s)]l 

P J 7t(s) J 


^ t 

Y - 7 w I/t'x' * w [ ds sin[ 7 '- 7 (s)]l 

* L M J M s) J 


Z - e 


■ 1(7*7* ) 


Substituting Eqs . (3.4)-(3.7) in Eq . (3.2), we obtain 


K(x,t; x',t) - ) 4 > (x'.t'^ 

[_ n n 


n-0 


where 


, . . .1 fiMTa 2 Z /•■' f* 1 . f(s) . . , .11 

n " eXP teU X • M ^ X J dS cos[7-7(s)]]J 

X eXp {ftf X ’ M j ds ^ 7 ^ sih[ 7 - 7 (s)]j 2 J 


x H n {^C 1 x ' 7~ [ ds sin[7-7(s]l| 

n 1 M L mv; J Ms) JJ 


x e 


i [ 0 ( t) - (n+ 2 > 7 ( t) ] 


Moreover, we may write 


V> n (x,t) - exp{i[(9(t) - ( n+“ ) 7 ( t ) ] ) ^ n (x,t) 


where 


^ n (x,t) 


_I_ M7 H h 
2 n n ! '* 




cos[7-7(s) ] 


]} 


(3.5) 

(3.6) 

(3.7) 

(3.8) 


(3.9) 

(3.10) 


356 



* ' « J d ‘ ^ sinl 


7-7(s)] 


“n - H I » s Sln[ ’- 1 ' (S > 1 }] ' 


(3.11) 


In Eq. (3.10), the wave function V> n (x,t) is merely a unitary transformation of 

<f> (x,t), and thus 4> (x,t) satisfies all the properties associated with V (x,t): 
n n n 


dx - <m I n> - dx </>*^> - 6 

J n n 1 J n n m,n 


(3.12) 


The expectation value of a given operator O is 


<m|o|n> - J dx V>* “ J dx ^n"^n 


(3.13) 


IV. Energy expectation values 

For the forced time - dependent harmonic oscillator system, both the Hamiltonian 
and Lagrangian have the units of energy but depend on time. We must therefore find 
a time - invar iant energy operator. If /3(t) is a particular solution of Eq . (2.3), 
we have 


(X-/3) + U 2 (t) (X-0) - 0 , 

dt 


and from Eqs. (2.26) and (2.27') we note that 


•»2 2 3 

rj + w ( t) rj - 0 /f? 


From Eqs. (4.1) and (4.2), we get the following expression for the energy; 


1 *2 ♦ •• M • • 2 M 2 x-B 2 

E - (r?p-M rjx) + (M«|X-«yp) + °<0i|-i|0) + J 0 (~ fi ) 


(4.3) 


Because Eq . (4.3) is time invariant, we can use it for the quantum mechanical 

energy operator, 

2 2 2 

„ M r? d M *2 • 2 . 2 nn. 2- 

E op - -m~l + 2^ + ” )x - 2 (2x a^ + X) 

r dx 


+ (Pb-nfa) + M^x) + Mf nV + | i,y 2 fi 2 + *(^-r//9) 2 


(4.4) 
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Equation (3.11) now simplifies to the expression 


2 2 


< (x.t) - f— +UX e' 45 (x ^> H 
'•2 n nt/jr 

- e^ +Bx Hf«(*-/J)] . 

^■2 n n ! Jn J n 


6 - (M7/J0 


(4.6) 


(4.7) 


" 1 


ds 1 1 cos[7-7(s)' 

7 7(s) 


0(t) - f ds sin[7-7(s); 

M/7 J 77 (s) 


A - ill - 6/2 


(4.8) 


(4.9) 


(4.10) 


B - iA + 06 


(4.11) 


Here, £(t) is a particular solution of Eq. (2.3). 

In order to evaluate the energy expectation E - <m|E |n>, we perform the 

oy ra , n op 

following calculations: 


x|n> ” y \s l n4 ’ 1> + 7n|n-l>] 

x 2 |n> - — (7(n+2) (n+1) |n+2> + (2n+l)|n> + /n(n-l) |n-2>] 


(4.12) 


(4.13) 


p | n> - ^ 2 (n+1) | n> + * B|n> + f(«+f) 72n|n-l> 

2 2 
p 2 |n> - -y 2 [^Y 7(n+2) (n+1) | n+2> +- 2/2 /H+I|n> + 2(A+^) (2n+l) |n> 

6 & 


+ 2/2 7^ + 5B)/n|n-l> + 2(^r+2A+5 2 )/n(n-l) |n>] 
5 5 2 
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iJ 


xp|n> - ^ ./(n+2) (n+1) |n+2> + y|j /n+T|n+l> + [^p2n+l)+n] | n> 

5 6 

+ y=J /n | n- 1> + (^+l)/n(n+l) |n-2>) (4.16) 

^ s 

px|n> - xp|n> + |n> , (4.17) 

Substituting Eqs. (4.6)-(4.17) into Eq. (4.4), we directly obtain the energy 
expectation values as 

- n)rf(rv^) . (4.18) 

This energy expectation value is a time- invariant quantity. 


V. Uncertainty values 

The uncertainty product defined as 

9 9 * o o U 

(AxAp)^ ^ - ([(<ra|x |n> - <m|x|n> ) (<m|x | n> - <m|x|n> )] 

x [(<m|p2|n> - <m|pjn>^) (<m|p^|n> - <m | p | n>^ ) ] ^ ) *^ . (5.1) 

Inserting Eqs. (4 . 12) - (4 . 15) into Eq . (5.1), we obtain 

; 2 4 i 

(AxAp) n n - (1 + 7 “-jjjr (n+pX 
' * ‘ 7 n 

;2 . 

(AxAp ) n+2 n - (1 + Tp-p 7 (n+2) (n+1) X 

’ 7 1 

2 . 

(AxAp) n n+2 - (1 + 7 ^ 7 ) Vn(n-l) X • 

’ 7 <1 

VI. Coherent states of the time-dependent harmonic oscillator 

First, we construct the creation operator a^ and destruction operator a. For 
a forced t irae- dependent harmonic oscillator, it Is not possible to construct a and 
a^ , but we can construct a and a^ for the time dependent harmonic oscillator. From 
Eqs. (4.12) and (4.14), we obtain 


(5.2) 

(5.3) 

(5.4) 
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(6.1) 


a t_ (fj) 4 [d + *7) x - i^] 
* m m 7 


a - (f/ td - 1^7) x + i^] 

f?7 M 7 


From Eqs . (6.1) and (6.2), we can represent (x.p) in terms of (a ,a) as 
x - (aUa) 


P - (“f)* 1 !<“: + i) a* * <a- - l).l . 

*n *n 


Also from Eqs. (6.1) and (6.2), if [x,p] - i\/i we see that 


[a 1 ,a] - 1 . 


Conversely, from Eqs. (6.3) and (6.4), if [a ,a] - 1 we note that (x,p) - i)A. 

The coherent state can be defined by the eigenstate of the nonherraitian 

_ 15 

operator a, 


a q> - a Q> 


Let us find the coordinate representation of the coherent state. From Eqs. (6.2) 
and (6.3), we have 


( 2w^ ( (1 " i h' ) x + *7 <x ‘ l 0> " Q<x ‘ l Q> • 

rj 7 m 7 

We solve this equation and change the variable x' into x for convenience, 


<x|a> - N exp + i^) x 2 + a x] 

n 7 ^ 


We choose the constant of integration N such that 


J dx | <x | a> | ^ - 1 


Then, we find the eigenvector of the operator a in the coordinate representation 
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( 6 . 10 ) 


, r Ma, 1 ,jLn 2 A /2M2.4 1| 1 2 1 2 

<x | q> - (^jj) exp[“(-l + i“) x + ( K ) ax ' 2 l a * " 2 a 


*?7 


Next, we show that a coherent state is a minimum uncertainty state. From Eqs. 

( 6 . 3 ), ( 6 . 4 ) and ( 6 . 6 ) and their adjoints, we evaluate the expectation values of x, 

2 2 

p, x and p In the state |ct>: 


<a|x|a> - ( a +a) 

2My 


( 6 . 11 ) 


< Q |p|a> - [(*7 + i)o* + - Do] 


V7 


*n 


( 6 . 12 ) 


<a|x Z |a> - “ ^ <a|a +a +aa +a a|a> 
2M7 


W *2 2 * 

- 1 **“ (a +a + 2 aa + 1 ) 

2M-y 


(6.13) 


<a|p 2 |a> - ^ {(*7 + i ) 2 a * 2 + (®7 • i ) 2 <* 2 


<77 


*77 


+ [ C 1 ^) 2 + 1 ] ( 2 aa*+l) ) 

n'y 


(6.14) 


The uncertainty value is 


2 2 * 

AxAp - [(<a|x |a> - <a | x | a>) (<a | p |a> - <a|p|a>)] 


- M/2 ll+C^) 2 ] 1 " 
*77 


(6.15) 


which is the minimum value allowed by Eq. (5.2), 


VII. Results and discussion 

In the previous sections, we have obtained the propagator, wave function, 
energy expectation values, uncertainty values and coherent state for a quantum 
forced time - dependent harmonic oscillator. These quantities represent the solution 
of the classical equation of motion for the time -dependent harmonic oscillator. If 
we set f(t) equal to zero, then our solution is correct for the time - dependent 
harmonic oscillator. Setting w(t) - u> q gives results for the forced harmonic 
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oscillator. For the case where f(t) - 0 and w(t) - w , our results are those of 

o 

the simple harmonic oscillator. 

For the explicit time - dependent system, we need to consider the quantum 

mechanical operator. In our work, the Hamiltonian, Lagrangian and mechanical 

energy have the units of energy, but these are not time Invariant. Yet, In order 

to solve macroscopic physical problems, we use time- invariant operators. For this 

reason, we have derived the energy operator from the classical equation of motion 

and used it to calculate energy expectation values. Our energy operator is similar 

1 2 

to the Ermakov-Lewis invariant operator. * Our quantum energy expectation vlaues 
are time - independent quantities, and our uncertainty values are consistent with 
Heisenberg's uncertainty principle. Yet, our uncertainty values are time 
dependent, in contrast to those of time - independent systems. 

Since it is not possible to construct a coherent state for the forced time- 
dependent harmonic oscillator, we have constructed It for the t ime - dependent 
harmonic oscillator. In general, the coherent state is a minimum uncertainty 
state, which is also true for our system. 

Time - dependent systems are observed in various physical experiments. Two 
general types of such systems are: that which is formed through its own 

environmental conditions, and that which is formed when external forces are added. 
In regard to the second type, various experiments are being carried out to see how 
an applied, time -dependent electric, magnetic of* other field can alter the physical 
properties of materials such as semiconductors and superconductors. Experiments 
show that a system becomes time dependent when a t Ime - dependent electric or 
magnetic field (such as a.c.) is applied. However, obtaining the quantum 
mechanical solution by a direct method is not easy mathematically. One way of 
obtaining a solution is to use the propagator method as indicated in this paper, 
where the relevant equations are those of a time - dependent harmonic oscillator. 

Our results, which are exact for one dimension, can be extended to two or more 
dimensions, and they can also be applied to time -dependent macroscopic systems. 
One example of an extension to two dimensions would be to solve the motion of a 
quantum electron in a time-dependent magnetic field. 
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